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Effect of the phase fluctuations of the order parameter on the stabiUty of the Fulde-Ferrell- 
Larkin-Ovchinnikov (FFLO) states are examined in exactly two-dimensional (2D) type-If super- 
conductors with cylindrically symmetric Fermi surface on the basis of a generalized Ginzburg- 
Landau theory. It is found that for the FFLO states with oscillations in a single direction, 
not only the long-range order but also quasi- long-range order (QLRO), which is characterized 
by a power law decay of the order parameter correlation function, is suppressed by the phase 
fluctuations at any finite temperatures. On the other hand, for the FFLO states with order 
parameter structures such as triangular and square lattices, it is shown that the QLRO is pos- 
sible as the uniform BCS state. Systems with anisotropy in the Fermi surface and pairing are 
also discussed. 
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tor, Kosterlitz-Thouless transition, phase fluctuations, Pauli paramagnetic limit, Chan- 
draseker and Clogston limit, quasi-long-range order 
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The nonuniform superconducting state called the 
Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state is char- 
acterized by spatial oscillations of the order parameter 
due to the Zeeman energy additional to the spatial de- 
pendence due to the vortices around the magnetic flux 
lines. It is known that the FFLO state is stable at high 
magnetic fields in type-II superconductors within a mean 
field theory if some ideal conditions are satisfied. "'^' ^ How- 
ever, there has not been any theory which has examined 
this state beyond the mean field approximation, although 
the two-dimensional (2D) systems are important when 
one considers the FFLO state, as we shall see below. In 
this paper, we examine the effect of the phase fiuctua- 
tions in the 2D systems on the stability of the FFLO 
state. 

For the FFLO state to occur in a type-II superconduc- 
tor, the following two conditions need to be satisfied. (1) 
Orbital pair breaking effect is sufficiently weakened by 
some mechanism, so that the Zeeman energy reaches a 
value of the order of the zero field gap, near the upper 
critical field. "^ (2) The sample is in a clean limit. ^' ^ In the 
alloy type-II superconductors, these conditions are diffi- 
cult to be satisfied, but it is possible to satisfy them in 
the organic and high-Tc oxide superconductors because of 
the strong electron correlations, narrow electron bands, 
and quasi-low-dimensionality. In low-dimensional super- 
conductors, the orbital pair breaking effect is weakened 
by applying the field in any direction parallel to the most 
conducting layer, ^"^^ with a sufficient accuracy. ^^ 

Further, the enhancement of the FFLO critical field 
due to the anisotropics of the Fermi surface (including 
the low-dimensional nature) and the pairing^^^" would 



increase the possibility of the FFLO state being found 
in the organic and oxide superconductors as discussed 
in detail in our previous papers. ^^'^^ This is not merely 
a quantitative issue, but makes a qualitative difference 
in the realization of the FFLO state, if one takes into 
account negative effects on the FFLO state, such as in- 
ternal field enhancement.^^^-'^^'^^ 

Therefore, low-dimensional exotic type-II supercon- 
ductors, such as the organic and copper oxide supercon- 
ductors, might be good candidates of the FFLO super- 
conductor. At present, we do not have any experimental 
evidence of the FFLO state, but some unusual experi- 
mental facts in organic superconductors can be explained 
as a consequence of the FFLO state, ^^^^^"^^ although 
other explanations such as the triplet pairing supercon- 
ductivity are also possible. ^^"^^ 

In this paper, we regard the quasi-one-dimensional 
organic superconductors such as (TMTSF)2X and 
(DMET)2X as quasi-2D systems, in the sense that (1) 
Fermi surfaces are sufficiently deformed by the applica- 
tion of pressure and so on, so that the nesting instabilities 
are suppressed, and (2) there are no flat portions in their 
Fermi surfaces, which would cause the FFLO critical fleld 
at T = to diverge. 

In many organic superconductors, especially in those 
with low zero-field transition temperatures to the su- 
perconductivity, interlayer electron transfers are strong 
enough to justify the mean field treatment at low tem- 
peratures of interest. In the theories of the FFLO state 
in the low-dimensional systems, ^^"'^^ such interlayer cou- 
plings are assumed implicitly. However, in some organic 
superconductors and the high-Tc oxide superconductors. 
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the interlayer coupling is not so strong. ^'^"'^^ For ex- 
ample, K-(BEDT-TTF)4Hg2.89Br8 has very small inter- 
layer coupling. The observed critical field reaches about 
5 times the Pauli paramagnetic limit (Chandraseker and 
Clogston limit) estimated from the zero-field transition 
temperature using a simplified formula.^'* The estimated 
value should be modified by many effects, but it seems 
difficult to explain such an extremely large value within 
the mean field theory if a singlet pairing is assumed. As a 
possibility, low-dimensional fluctuations may play an im- 
portant role in the large ratio of the upper critical field 
to the zero field transition temperature. 

In such low-dimensional superconductors, it is known 
that the phase fluctuations of long wavelength would 
suppress the long-range order (LRO) of the supercon- 
ductivity at any finite temperature, but the quasi-long- 
range order (QLRO) may occur below a finite Kosterlitz- 
Thouless (KT) transition temperature.'^^' "^^ Here, the 
QLRO is characterized by a power law decay of the order 
parameter correlation function for long distances. When 
small interlayer transfer is introduced, the LRO of the su- 
perconductivity occurs at the temperatures below a finite 
transition temperature close to the KT transition tem- 
perature, which would be much smaller than the mean 
field value of the transition temperature. Also for the 
FFLO state, the phase fluctuation must have a serious 
effect on the LRO in the low-dimensional systems. In 
two dimensions, it is clear that the LRO of the FFLO 
superconductivity does not occur at finite temperatures, 
because of the existence of continuous symmetry, as in 
the BCS superconductivity. However, it is not yet clear 
whether the QLRO of the FFLO superconductivity oc- 
curs, since the nonuniformity of the FFLO state may 
aff'ect the phase fluctuations. 

We first note that the FFLO state could take various 
structures other than simple oscillations in a single di- 
rection at high magnetic fields. "'^^ The most general form 
of the order parameter ■)/;(r) is written as 



V-Cr) 



J2^t^e^^-- 



(1) 



near the FFLO critical field, where q^'s are some of the 
optimum wave vectors which yeild the highest critical 
field. For the cylindrically symmetric Fermi surface, we 
have an infinite number of optimum vectors for s-wave 
pairing, while we have four such vectors for d-wave pair- 
ing. We write their magnitude as g = |qm|, because they 
are the same. In this paper, since we consider only peri- 
odic solutions, only a finite number of qm's are consid- 
ered in eq. (1). The FFLO state with spatial oscillations 
in a single direction is expressed by eq. (1) in which only 
a single wave vector qi is considered, or two antiparal- 
lel wave vectors (qi = — q2) are taken into account. We 
call such structures one-dimensional (ID) structures. On 
the other hand, when we consider more than one q^ 
in eq. (1) which includes two q^'s neither parallel nor 
antiparallel to each other, we call it a 2D structure. In 



the FFLO state, the spatial symmetry is broken sponta- 
neously on the appearance of the order parameter. 

In the mean field approximation, the optimum form 
among those expressed by eq. (1) is determined by the 
terms of the fourth order of V'(r) in the free energy.^' ^^ 
It has usually been believed that a ID structure is opti- 
mum, since Larkin and Ovchinnikov found that it is opti- 
mum in a spherically symmetric system,^ but in practice, 
it may not be the case in other systems. For example, in 
a cylindrically symmetric system, the triangular, square, 
and hexagonal states occur as the temperature decreases 
for the s-wave pairing, while the square states occur at 
low temperatures for the d-wave pairing. ^^ This is not 
merely due to a geometrical effect, but due to the high 
FFLO critical field in low-dimensions, for which the 2D 
structures are favored over the ID structures due to spin- 
polarization energy. 

In general, the optimum structure may depend on the 
symmetries and other details of the Fermi surface and the 
pairing, and the strength of the field. Thus, we consider 
the general form of eq. (1) as the stable structure of the 
order parameter in the absence of the phase fluctuations, 
assuming that it is known a priori. 

We first construct a phase Hamiltonian from a gener- 
alized form of the Ginzburg-Landau (GL) Hamiltonian. 
We assume the rotational symmetry of the system and 
define the Hamiltonian as 



H[i^] 



dV [icIWWP 



1 



d|AV'(r)|^ 



(2) 



where we have taken into consideration only second or- 
der terms related to the phase fluctuations. The sec- 
ond term has been added so that the energy is bounded 
when c < 0. This is the most compact form which al- 
lows nonuniformity of the order parameter in the ab- 
sence of the vector potential. The detailed form of the 
GL free energy near the tri-critical point was microscopi- 
cally calculated by Buzdin and KachKachi,^^ but we take 
a more phenomenological viewpoint. Substituting eq. (1) 
into eq. (2), we obtain an expression for the energy as a 
function of q. When c < 0, a finite value oi q = y^\c\/{2d) 
gives the energy minimum, and the nonuniform state is 
stabilized. 

For finite temperatures, we introduce the phase fluc- 
tuation in the order parameter '>p{r) by 



ip{r) = iP{r)&'^ 



(3) 



where V'(r) is the order parameter at T = given by 
eq. (1) with the above value of q, and (j)(r) is a real func- 
tion expressing the phase fluctuations. We neglect fluc- 
tuations in the magnitude of the order parameter at suf- 
ficiently low temperatures. Assuming that (/)(r) is small 
and varies slowly over a distance scale much larger than 
the inverse of q, the phase Hamiltonian is derived as 



^0 = ^rfEl^o"^l'/dM(A0)^+4(q„ 



V4>n (4) 
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When all values of tp^ are equal to ipo, 

H^ = C^ ld\[{^4>f + ±^(q„ . V0)2], (5) 

where Nq denotes the number of ?n's which are taken 
in the summation in eq. (1), and C^ = ^d ■ Nq\tpo\^. 
Here, we find a remarkable difference in the second term 
between the f D and 2D structures: for the f D structures, 
the second term is proportional to (q • V(/))^, while for 
the 2D structures, it is proportional to g^(V(/))^. 

We should note that the FFLO states have many line 
nodes in the real space, which divide the real space into 
many sections or cells. Since the phase (j) is not actually 
defined where the amplitude vanishes, it appears that the 
phase can be discontinuous there. However, if the phase 
jumps between the sections divided by the line nodes, 
Josephson current is induced and causes the phase dif- 
ference to vanish rapidly. If we regard the connection 
between the sections as a tunnel junction, it is almost 
transparent to electron tunnelling, because no barrier ac- 
tually exists between the sections. Thus, the phase differ- 
ence must be negligibly small even if it exists. Therefore, 
it is reasonable to assume that the phase function (/) is a 
smooth continuous function over the sections. 

The spatial correlation of the order parameter is de- 
fined by 

C(r) =<[^(r)]*^(0)> (6) 

with tj){r) = ^o(i^^ ■ Since eq. (5) is bilinear in (/>(r), we 
have 



C(r) 



|^opexph-<[0(r)-(/)(O)]2>] 



lV'opcxp[ 



-ij5:|e*-r-iP<|<^,P>], 

k 



(7) 
where L is the linear dimension of the system. The aver- 
age <|(/'kP> is calculated from eq. (5) as 



< 



'>= 



Cd, 






(qm • k) 



(8) 



The states with the ID structures are given by A''^ = I 
and qi = q, or A^^ = 2 and qi = — q2 = q. For the both 
states, we find 



<\(t>u 



->= 



I 



C^ k^ + - 



(9) 



4(q.kf 

We substitute this expression into the integral in eq. (7) . 
Unless k _L q, the integrand is proportional to l/k'^ for 
small fc, and then the integral over k diverges logarithmi- 
cally. However, in practice, we have a stronger divergence 
from the contribution near k _L q. Since the essentially 
important contribution comes from the vicinity of k _L q, 
we put (f = 7r/2 in the factor je*'"" — Ip, where (p is the 
angle between k and q, except when the angle 9 between 



q and r is not very small. Then, this factor introduces 
the cutoff momentum of the order of I/(r sin 9) for small 
k. For long distances such as r ^ ^ ~ l/q, the factor 
jgikr _ i|2 (.^jj^ j^g replaced by the average value of 2. 
Thus, we obtain the short-range order 

C(r)-|V;opexp(-r/e) 



(10) 



with ^ = 2TTqCcf,/{Tsm9). When r//q, since the in- 
tegrand vanishes at (p = 7r/2 because of the factor 
l^ik-r _ ^2^ .^p obtain a logarithmic divergence of the 
integral in the long distance limit and the power law de- 
cay of the correlation function. Except for this direction 
r/^q, we only have a short-range order for the FFLO 
states with the ID structures. 

On the other hand, for the states with 2D structures 
such as the triangular lattices and the square lattice, ^^ 
a QLRO occurs in two dimensions. In the denominator 
of the integrand of eq. (8), the terms in the summation 
Sm(l"i ' ^)^ '^'^ ^"^t simultaneously vanish when k ^ 0. 
Hence, for long wavelengths, the first term in the denom- 
inator proportional to k* can be omitted, and we have 
a logarithmic divergence of the integral in eq. (7) , which 
leads to the power law decay of the correlation function in 
the long distance limit. There are three possible 2D struc- 
tures, that is, the triangular state: Ng = 3, qi = (j(I,0), 
q2 = g(-I/2,V3/2), qa = g(-I/2, -V3/2), the square 
state: Nq = 4, qi = (7(1,0), q2 = (7(0,1), qs = -qi, 
q4 = — q2, and the hexagonal state: Nq = 6, qi = q{l, 0), 
q2 = q{-l/2,V^/2), qg = g(-I/2, -V3/2), q4 = -qi, 
qs = — q2, qe = — qa- For these states the Hamiltonian 
is reduced to 

H^ = ]^K^ jd?r{V<pf (II) 

with K(j, — 2q^C(f,. The summation ^j^ in eq. (7) is taken 
between cutoff momenta k ^ \/r and k ~ 1/a, where a is 
a length of the order of the lattice constant. The former 
cutoff is due to the factor je'''''— Ip. Therefore, we obtain 
a power law decay of the correlation function 



C(r) = |^o|'(-)'' 



(12) 



with f] = T/{2ttK^) = T/{2TT\c\Nq\tPo\^). The KT cri- 
terion for the stability of the QLRO against free vortex 
formation is roughly written as 

T<^K^ = ^\c\Nq\ijo\^^TKT. (13) 

Since the stiffness constant K^f, is reduced by the polar- 
ization of vortex pairs, the actual value of Tkt is smaller 
than that given above. In the magnetic field, the vor- 
tex formation would reduce the energy by an amount 
proportional to the number of the vortices, since they 
are accompanied by the spin and orbital polarizations. 
Thus, the KT transition temperature decreases with the 
magnetic field, and vanishes at a critical magnetic field. 
The phase diagram will be investigated in a future study. 
We briefly comment on the phase fluctuations in three 
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dimensions. In this case, the spatial correlation of the 
phase diverges logarithmically for the FFLO states with 
ID structures. Thus, it may be conjectured that the 
LRO's of the FFLO states with ID structures are sup- 
pressed by the phase fluctuations at finite temperatures 
even in three dimensions, and a QLRO might occur. This 
subject will be discussed in more detail in another paper. 
Lastly, we discuss anisotropic systems. We have stud- 
ied the isotropic Hamiltonian given by eq.(2) in two di- 
mensions, but it would be interesting to extend our cal- 
culation to anisotropic systems for more realistic models. 
For example, a model Hamiltonian 



Hm = Jd\ [^c|v^(i 



)P + i(cio + 3di)|A^(r)|^ 



"diEJ(e.-V)V(r)P] 

(14) 
describes a system with tetragonal symmetry, where ei = 
-64 = (1,0), 62 = -63 = (0,1) and do > ^i > 0. For 
an order parameter of the form "tpir) = -00 exp(iq • r), the 
energy of the system is given by 

E{q, ip) = L^li^ol^i^cq^ + ^q^{do - d, C0s(4^)}], (15) 

which becomes minimum at ip ~ 0, 7r/2, tt, 37r/2 and 



q = ^\c\/2ido-di). 
The phase Hamiltonian is derived as 

+2{do - 3di)(q,„ • V0)2 + i(do + 3di)(A0)2 

-2d^m^r + id^cj^n] 

(16) 
It is easily verified that the QLRO's of the FFLO states 
with ID structures are stabilized as a local free energy 
minimum, differently from the isotropic model, but the 
stiffness constants are smaller than that in the state with 
a 2D structure, i.e., a square state in this case.^^ There- 
fore, the KT transition temperature to the FFLO state 
with the ID structures is smaller than that to the state 
with the 2D structure also in anisotropic systems. 

The results of this paper are explained as follows. In 
ID structures, the phase fluctuations with small wave 
vectors k perpendicular to q are regarded as the fluc- 
tuations of the direction of q. On the other hand, in 
2D structures, any phase fluctuations are inevitably ac- 
companied by fluctuations of the magnitude of q^'s in 
addition to the fluctuations of the directions. Thus, the 
stiffness constants in the 2D structures are usually larger 
than those in the ID structures, and the phase fluctua- 
tions affect the FFLO states with the ID structures more 
seriously. 

In conclusion, we have examined the effect of the phase 
fluctuations on the FFLO state in exactly 2D type-II su- 
perconductors. From the generalized GL Hamiltonian in 



which a nonuniform state is stabilized at T = even 
in the absence of the vector potential, we have derived 
an effective Hamiltonian for the slowly varying phase of 
the order parameter. On the basis of the phase Hamilto- 
nian, we have obtained the following results. (1) When 
the FFLO state at T = has an order parameter which 
oscillates in a single direction, such as A ^ cos(q • r) 
and A ^ exp(iq • r), the LRO and even the QLRO's are 
unstable at finite temperatures because of the phase fluc- 
tuations in the directions perpendicular to q. (2) When 
the FFLO state at T = has an order parameter with a 
2D structure such as the triangular and square lattices, 
since the phase Hamiltonian has the same form as in 
the case of uniform superconductivity, a QLRO occurs 
below a finite KT transition temperature, although the 
LRO is suppressed. Therefore, in 2D FFLO superconduc- 
tors with very weak interlayer coupling and cylindrically 
symmetric Fermi surface, the structure of the order pa- 
rameter must be that of the 2D lattices, and not simple 
oscillations in a single direction. 
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